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1 Introduction 



Most classical statistical mechanics [1, 2] and two-dimensional quantum field theory models [3] exact 
solubility relies on the 7?.-matrix TZ{u), where w is a spectral parameter, acting on the tensor product 
V of a given vector space V . lZ{u) abridges the Boltzman weights of classical statistical mechanics 
models or the two-body scattering matrix of integrable quantum field theories and is a solution of the 
Yang-Baxter (YB) equation 

7^12(«)7^13(M + v)nM = 7^23(v)7^l3(« + v)7^i2(«), (i.i) 

in V^®V'^®V^ , where TZi2 = TZ^l, 7^23 = f^'R-, etc. 7^(u) may also depend on an additional parameter 
r] in such a way that 



n{u,7)) = l + r]r{u)+o{'rf), (1.2) 

where 1 is the identity operator in the space V ®V . In this case the solution 72.(u) of Eq.(l.l) is said to 
be semiclassical and the "classical r-matrix" obeys the equation 

[ri2{u),rx2,{u + v) + r23{v)] + [r\z{u + v),r23{v)] = 0. (1.3) 

called the classical Yang-Baxter equation. Nondegenerate solutions of (1.3) in the tensor product of two 
copies of a simple Lie algebra g , rij{u) G gi ®gj , i,j — I, 2, 3, were classified by Belavin and Drinfeld [4] 
and its important role in the theory of classical completely integrable systems can be found, for instance, 
in [5]. 

In the skew-symmetric case rji{—u)+rij{u) = 0, the classical YB equation embodies the compatibility 
condition for the system of linear differential equations 

=Y, nj{z, - Zj)^{z„ ...,^;v) (1.4) 

in TV complex variables zi,...,zn for a vector- valued function ^'(zi, zjv) with values in the tensor 
product space V = V'^ (E) ■ ■ ■ (S) . k is a yet non specified coupling constant. 

In the rational case [4]. very simple skew-symmetric solutions are known: r{u) = C2/U, where C2 G 
g (8" g is a symmetric invariant tensor of a finite dimensional Lie algebra g acting on a representation space 
V, and the above system of linear differential equations (1.4) is known as Knizhnik-Zamolodchikov (KZ) 
system of equations for the AT— point conformal blocks of WZW conformal field theories on the sphere [6] . 
The trigonometric solutions of Eq.(l.l), or quantum group solutions, were classified by Jimbo [7], and in 
this case the system of equations (1.4) is named generalized system of KZ differential equations. 
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The Gaudin [8] Hamiltonians Gi are related to the classical r-matrices by 

Gi = ^rij{zi - Zj) (1.5) 

and the condition for their commutativity is again the classical YB equation (1.3). 

This interplay of the classical YB equation with conformal field theory and the theory of Gaudin 
Hamiltonians can be understood through the algebraic Bethe Ansatz, here shortly described in the 
following way. 

Babujian and Flume [9] showed that the algebraic Bcthc Ansatz [10] for the theory of the Gaudin 
models [8] can be gotten from the the algebraic Bcthc Ansatz of a non-homogeneous lattice vertex model, 
the appropriately expanded body Bcthc wave vectors rendering, in the scmiclassical limit, solutions 
to the generalized system of KZ equations. For instance, in the rational su{2) example, the algebraic 
quantum inverse scattering method [10] allows one to write the following equation 

V— -\ J' *1**^ 

t{u\z)^{ui^... Up) = A{u,ui, ■ ■ ■ ,Up\z)^{ui, ■■■ , Wp) - ^ — ^ • (1.6) 

a=l 

Here t{u\z) denotes the transfer matrix of the rational vertex model in an inhomogeneous lattice acting on 
an TV-fold tensor product of su{2) representation spaces. meaning = <I>(ui, • • • Ua-i,u, Ua+i, Up) 
; J-'a{ui,--- ,Up\z) and A{u,ui, ■ ■ ■ ,Up\z) being c number fmictions. The vanishing of the so-called 
unwanted terms, J-a = 0, is enforced in the usual procedure of the algebraic Bethe Ansatz by finding the 
parameters ui, ...,Up. In this case the wave vector "I>(ui, • • • ,Up) becomes an eigenvector of the transfer 
matrix with eigenvalue A(u,wi, • • • ,Up\z). If we keep all imwanted terms, i.e. J-'a ^ 0, then the wave 
vector $ in general satisfies the equation (1.6), named in [11] as off-shell Bethe Ansatz equation (OSBAE). 

There is a neat relationship between the wave vector satisfying the OSBAE (1.6) and the vector-valued 
solutions of the KZ equation (1.4): the general vector valued solution of the KZ equation for an arbitrary 
simple Lie algebra was found by Schechtman and Varchenko [12]. It can be represented as a multiple 
contour integral 



\E'(^i, . . . , ^jv) ~ J ' " J '^('"1' ■■■,Up\z)(p{ui, ...,Up\z)dui ■ ■ ■ dup. (1.7) 

where the complex variables z\,...,zn of (1.7) are related with the disorder parameters of the OSBAE . 
The vector valued function <f){ui, ...,Up\z) is the scmiclassical limit of the wave vector $(ui, ...,Up\z). In 
fact, it is the Bethe wave vector for Gaudin magnets [8], but "off -shell ". The Bethe Ansatz for the 
Gaudin model was derived for any simple Lie algebra by Reshetikhin and Varchenko [13]. The scalar 
function X{ui, ...,Up\z) is constructed from the semiclassical limit of the A('u = Zk\u\, ...,Up\z) and the 
•?^a(wl, • • • ,Up\z) functions. 
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In this paper we investigate the semiclassical hmit of Izegin-Korepin (IK) model [14], which corre- 



(2) 

spends to the twisted affine Lie algebra A2 solution of [7]. 



(2) 



The paper is organized as follows. In Section 2 we present the algebraic Bethe Ansatz for the A2 
vertex model. Here the inhomogeneous Bethe Ansatz is read from the homogeneous case previously 
derived for the 19- vertex models [15, 16]. We also derive the off-shell Bethe Ansatz equation for this 
vertex model. In Section 3 , taking into account the semiclassical limit of the results presented in the 
Section 2, we describe the algebraic structure of the corresponding Gaudin model. In Section 4, data of the 
off-shell Gaudin equation are used to construct solutions of the trigonometric KZ equation. Conclusions 
are reserved for Section 5. 



2 Inhomogeneous Algebraic Bethe Ansatz 

Besides 7?. , it is usual to consider the matrix R = VTZ that satisfies the equation 

i?12(w)i?23(w + v)Ri2iv) = i?23 (v)i?12 (w + v)R23{u). 



(2.1) 



We choose the regular solution of the YB equation for the IK 19-vertex model [14] to be normalized 
such that 
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(2.2) 



where 



Xi (u) = sinh(u -|- 2r}) cosh(u -|- ^rf) , X2 (w) = sinh u cosh(M + 3»7) , 

xs{u) = sinhucosh(M + 77), X4{u) = sinhMC0sh(M + 3?7) + sinh 27? cosh Sr/, 

X5{u) = e" sinh2?7Cosh(w + 3?7), y5{u) = c^^ sinh2r] cosh{u + Si]), 

X6{u) = e""''^'' sinh 277 sinh u, yeiu) = —e~"~^'' sinh 2?7 sinh m, 

x^{u) = e2"xi(«)-e2»+4''x3(u), yriu) = e-^^'xiiu) - e'^^'-^^'xaiu) 



(2.3) 



An inhomogeneous vertex model has two parameters: a global spectral parameter u and a disorder 
parameter z, so that the vertex weight matrix TZ depends on the difference u — z and consequently the 
monodromy matrix defined bellow will be a function of the disorder parameters Zi. 



3 



The quantum inverse scattering method is characterized by the monodromy matrix T(u\z) satisfying 
the equation 



R{u - v) [T{u\z) (8) T{v\z)] = [T{v\z) O T{u\z)] R{u - v), 



(2.4) 



whose existence is guaranteed by the YB equation (1.1). T{u\z) is a matrix in the space V (usuahy called 
auxiliary space) whose matrix elements axe operators on the states of the quantum system (which will 
also be the space V in this work ). The monodromy operator T{u\z) is defined as an ordered product of 
local operators (Lax operator), on all sites of the lattice: 



T{u\z) = Cn{u - zn)jCn-i{u - zn-i) ■ ■ ■ Ci{u - zi). 
We normalize the Lax operator on the n*'' quantum space to be given by: 



(2.5) 
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-^22 {u — Zn) ^23 i^^ ^ Zn) 
L33 - Zn) 



(2.6) 



L^^^{u), a,j3= 1, 2, 3 are 3 by 3 matrices acting on the n*'* site of the lattice, so the monodromy matrix 
has the form 



Ai(m|z) Bi{u\z) B2{u\z) 
T{u\z)=\ Ci{u\z) A2{u\z) B^{u\z) 
C2{u\z) CMz) AMz) 



(2.7) 



where 



TiMz) = 



ki 

1,2,3. 



3 

E 

.,fejv-i=l 



' 4ffe2^^(" - ^Jv-i) (8> • • • (8) L^kLiA'^ ~ ^i)- 



(2.J 



The vector in the quantum space of the monodromy matrix T{u\z) that is annihilated by the operators 
Tij{u\z), i > j {Ci{u\z) operators, i = 1,2,3) and it is also an eigenvector for the operators Tii{u\z) ( 
Ai{u\z) operators, i = 1,2, 3) is called a highest weight vector of the monodromy matrix T{u\z). 
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The transfer matrix t{u\z) of the corresponding integrable spin model is given by the trace of the 
monodromy matrix in the space V 

t{u\z) = Ai{u\z) + A2{u\z) + A3{u\z). (2.9) 

The algebraic Bethe Ansatz solution for the inhomogeneous IK vertex model can be obtained from 
the homogeneous case [15]. The proper modification is a local shift of the spectral parameter u ^ u — Zi, 
and the functions needed in the sequel are 



X2{u) sinhu ' 2/6 (w) sinh277 ' 

, , Xi{u)x3{u) cosh.{u + rj) 

Xi{u)x3{u) — XQ{u)ye{u) cosh(u — ??) ' 

We start defining the highest weight vector of the monodromy matrix T{u\z) in a lattice of N sites as 
the completely unoccupied state 

|0>=®f=i( I . (2.11) 

Using (2.8) we can compute the normalized action of the monodromy matrix entries on this state 

Ai{u\z)\0) = X,(u|z)|0), C\{u\z) \0)=0, Bi{u\z)\0)^{0,\0)}, 

X,{u\z) = n ""/' V i = 1,2,3. (2.12) 

The Bethe vectors are defined as normal ordered states 4'„(wi, • • • , m„) which can be written with aid 
of the recurrence formula [15]: 

^'„(«i, ...,Un\z) = Bi(mi]^;)*„_i(u2, ■■■,Un\z) 

-B2(ni|z)f^^i%^ f[ Z{Uk-Uj)<i>r,.2{U2,:;Uj,...,Un\z), (2.13) 

with the initial condition = |0) , 'J'i(miI^) = Bi{ui\z) |0). Here Uj denotes that the parameter Uj 
is absent: ^{uj\z) = '^{ui, ...,Uj-\,Uj+i, - ■ ■ ,Un\z). 

The action of the transfer matrix t{u\z) on the Bethe vectors gives us the following off-shell Bethe 
Ansatz equation for the A2 vertex model 
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n n J — 1 

j=i j=2 1=1 (2.14) 

We now briefly describe each term which appear in the right hand side of (2.14) ( for more details 
the reader can see [16]): in the first term the Bethe vectors (2.13) are multiplied by c-number functions 
A„ = A„(u,ui, ...,Un\z) given by 

The second term is a sum of new vectors 
multiplied by c— number functions J-'j^ given by 

n n 

= X,{uj\z) H Z{uk - uj) - X2{uj\z) H Z{uj - Uk). (2.17) 
Finally, the last term is a coupled sum of a third type of vector-valued functions 

*Jl-2)=^2(«|^)M/n-2(«i,W,), (2.18) 

with c— number coefficients 

n 

^//"^^ = GijXi{ui\z)Xi{uj\z) Z{uk - ui)Z{uk - Uj) 

n 

+YijXi{ui\z)X2{uj\z) Yi ^i'^k - ui)Z{uj - Uk) 

k=l,kjtj,l 

n 

+FijXi{uj\z)X2{ui\z) Yl Z{ui - Uk)Z{uk - Uj) 
k=i,kjtj,l 

n 

+HijX2{ui\z)X2{uj\z) H Z{uj-Uk)Z{ui~Uk). (2.19) 
k=i,kjtj,i 
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where Gij , Yij , Fij and Hij are additional functions defined by 

x-jiui — u) 1 z{ui — u) xc^{uj — u) 1 



Gij = 



X3{ui - u) y{ui - Uj) oj{ui - u) X2{uj - u) y{u - ui) ' 
y'j{u-ui) 1 y^{u-ui) 1 



xs{u - ui) y{ui - Uj) X3{u - ui) y{u - Uj) ' 
V 1 / / \y5{u-Uj) ysiu-ui) y5{ui-Uj)\ 

y{u-Ul) [ X2{u-Uj) X2{U - Ul) X2{UI - Uj) ) 



^ Vbju - Ul) \ y^{ui - Uj) 1 ^ ^(m-m;) 



X2{u - Ul) \ X2{ui - Uj) y{u - Ul) w{u - Ul) y{u - Uj) 

_x_^{u-ui)_ 1 \ (2.20) 

X2{u - Ul) y{ui - Uj) J 

In the usual Bethe Ansatz method, the next step is to be impose the vanishing of the so-called unwanted 
terms of (2.14) in order to get an eigenvalue problem for the transfer matrix. 

We impose — and = into (2.14) to recover the eigenvalue problem. This means that 

4'„('Ui, ...,Un\z) is an cigcnstatc of t{u\z) with eigenvalue A„ , provided the parameters Uj are solutions 
of the inhomogeneous Bethe Ansatz equations 

N n , _ . 

l[z{Uj-Za) = n !(2-!V ^"'~"^"^' 

a=l fc=l, k^j 

j = l,2,...,n. (2.21) 
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3 Structure of the A^^^ Gaudin Model 

In this section we will consider the theory of the Gaudin models and to do that we need to compute the 
expansions of some results presented in the previous section in a power series of . 

As we already mentioned the IK vertex model corresponds , in the Jimbo [7] quantum group classifi- 
cation, to the twisted afhne Lie algebra A^^ ' ,which is constructed out of the second order automorphism 
of the Lie algebra , the swapping of its simple roots, such that the fixed algebra under this automor- 
phism, say go, is isomorphic to the Lie algebra Ai. We recall [17] the following base of "twisted " su{3) 
generators in the fundamental representation 

H = 



S+ 



W+ = \ -1 \ , W- = \ 1 (3.1) 



H, t/+ and U correspond to the oven sector of the automorphism and generate qq , the remaining 
generators correspond to the odd sector and transform according to the "isospin " 2 representation of 
50- 

The quadratic Casimir operator is 

C2=H^ + ^Y^ + 2{S+, S-} + {U+, U-} + {W+, W-}. (3.2) 

where {A, B} = AB + BA. 

In order to expand the matrix elements of T{u\z), up to an appropriate order in rj, we will start by 
expanding the Lax operator entries defined in (2.6): 

smh(2u — 22;„) \2 4 cosh(u — 2;„)^/ 

smh(2M — 2znj \4 cosh(w - Zn) ) 

- 1 I , J(2 + yn)-^nCOsh(2.-2.„) fx 3 g^ + g„ \ 

33 - ^ + ^'' sinh(2u-2^„) 1^2''" 4cosh(u-^„)2j ^• 





(3.3) 
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and for the elements out of the diagonal we have 

Lin) _ W-+e^("--)L/- , j-(n) W+ + e-^("-^-^)U+ 



sinh(2u — 22„) sinh(2u — 2z„j 



= 27? 0(7,2), =27? . ,,5" , . +o(r;2). (3.4) 

' sinh(2u - 22„) ^ ' ^' ' sinh(2u - 22„) ^'^ ^ ^ 

Substituting (3.3) and (3.4) into (2.8) we will get the semiclassical expansion for the monodromy 
matrix entries. For the diagonal entries we get 

A,{u\z)^l + 27jA[^\u\z) + Ar]^A\^\u\z) + o{r]^), i ^ 1,2,3. (3.5) 

where, the first order terms are given by 

and the second order terms are 

^ 1 Q 172 _ IT 

= E( 2 ^° ~ 4 cosh(L - 2I )2 ) + E "°^^(27x - 2z,) coth(27x - 2z,)g, ^ H, 

a=l ^ ""^ a<6 

^ ^ 1(2 + Yg) 1(2 + n) ^ ^ 1(2 + y„) gbC0sh(2M - 2^fc) 
^ sinh(2u - 2za) sinh(2u - 2zb) ^ sinh(2u - 2za) sinh(2u - 22:6) 

cosh(27i-2z„)ga(8)|(2 + n) ^ ^_ S- (E) S+ 



sinh(2u — 2za) sinh(2u — 2zh) sinh(2'u — 22;a) sinh(2w — 2zi,) 

sinh(2u — 2za) sinh(2tt — 2zb) ' 

a<b 



A^s\u\z) = f2ilH',-^-^^^^) + J2<^othi2u-2z,)coth{2u-2zb)H,^Hb 



a— 1 ^ ^ a<o 



^ ^ 1(2 + Yg) ^(2 + n) ^ 1(2 + y,) ® Ht,cosh{2u - 2zb) 

sinh(2'U — 2za) sinh(2u — 2zb) sinh(2u — 2za) sinh(2u — 22:6) 

^ cosh(2M - 2za)Ha g(2 + n) sr^4^_ S+^Sj; 



-E 



sinh(2u — 2zo) sinh(2u — 2zb) sinh(2u — 2za) sinh(2u — 2zh) 

-+ _ e-2(„-..)^+ x ^/ Wj- e2("-^^)t/, 
sinh(2u — 2za) J \ sinh(2u — 2zb) 
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JV 



3 Ya-Hl 



^_ |(i-i;)cg,|(i-rb) 

^—^ sir 



\4 cosh(2u — 2za)'^ ) sinh(2u — 2za) sinli(2M — 2zb) 

sir 



a<b 



sinh(2u — 2za) sinh(2u — 2zb) 



i!-^ sir 



a<6 



sinh(2u — 2za) sinh(2u — 2zb) 



It will be only needed to expand the off-diagonal elements up to the first order in r] 

Biiu\. 



B2{U 



sinh(2u — 2za) 



N 



2S: 



sinh(2u — 2za) 

;inh(2.-2..) )' 



' SI 



sinh(2u — 2za) 



- + o(r?^), 



N 



' .il^ SIT 



2St 



sinh(u - Za) 

a=l 



sinh(2u — 2za) 

a=l 

^Prom these expansions we have the following expansion for the transfer matrix (2.9): 



AT 



N 



t{u\z) = 3 + 277^ 



JV 



sinh(2u — 2za^ 
2 



2 cosh(u — Za)"^ ' 



a<b 



sinh(2u — 2za) sinh(2u — 2zb) { 3 



- + Ha '§> Hb cosh(2u — 2za) cosh(2'u 



+ ^Ya ® n + 2{S+ ® + S- S+) + W+ ® W,- + W- w+ 
+ e2^a-2^oc/+ (g) [/- + e-2^<.+2^!'[/- (g) ^7^} 
= 3 + 2r]T'-^\u\z) + 4:jfT'-'^\u\z) + o{rf). 
The second order contribution for t{u\z) is 



N 



N 



T^''\u\z)=Y.^ab{u)+Y.^Hl 
a<b 



which, with aid of the identity 

1 



a=l 



1 



2 cosh(u — Za)^ ' 



(- 



-2(ti-z„) 



-2(tt-2zi) 



sinh(2u — 2za) sinh(2u — 2zb) sinh(2za — 2zb) sinh(2-u — 2za) sinh(2u — 2zb) 



:) 
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can be written in the form 



where 



T^^\U\Z) - J2 g2(„_.J ^^1(1^ _ 2 ) + E(^a 2 COSh(l. - Za)^ ^ ^^'^^^ 
a=l ^ ' o=l ^ ' 



1 f 2 

Ga (u) = V ^-tt:^ TT^ { :j + cosh(2u - 2z„) cosh(2w - 2zb)H^ ® 

fr" smh(22a - 2zb) y 3 



+ \y^®Y, + 2 {St 5,- + S- ® 5+) + W^- W+ + W+ 



b 



(3.13) 



Here we observe that Qa{u) is nothing but the sum of semiclassical trigonometric r-matrices. This fact 
follows from the construction of the semiclassical r-matrices out of the quadratic Casimir: 

C2 = + 1^2 + 2{S-S+ + S+S~) + {U-U+ + U+U-) + {W+W- + W-W+). 

^ (3.14) 

The Gaudin Hamiltonians are defined as the residue of t{u\z) at the point u = Za- This results in N 
non-local Hamiltonians 

N 



1 r 2 

^« = g sinh(2.„ - 2z,) [ 3 + - ^ 



+^Fa ® n + 2 {St ® S^ + S- ® St) + W- w+ + wt ® w^ 
+ e2(^<'-^'>)i7+ ® C/- + e-2(^<»-^'')?7- C/b^} , 
a = 1,2, ...,7V. (3.15) 



satisfying 



In the next section we will use the data of the Bethe Ansatz presented in the previous section in order 
to find the exact spectrum and eigenvectors for each of these — 1 independents Hamiltonians. 

We complete this section deriving the A2 Gaudin algebra from the semiclassical limit of the fun- 
damental commutation relation (2.4): the semiclassical expansions of T and R can be written in the 
following form 

T{u\z) = 1 + 2r][l{u\z) + (3{u\z)] + o{r]'^), R{u) = P[l + 2r][r{u) + /?(«)] + 0(7/^)] . 

(3.17) 
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where 



JV 



3 sinh(2'u - 2za) 

a=l ^ ' 

Using (3.7-3.8) one can see that the "classical Z-operator " has the form 

/ y{u\z) + n{u\z) W-{u\z)+U-{u\z) S-{u\z) 



l{u\z) 



\ 



W+{u\z) + U+{u\z) -2y{u\z) -W^{u\z)+U-{u\z) 

s+{u\z) -yv+{u\z)+u+{u\z) y{u\z)-n{u\z) J 



where 



N 



N 



n{u\z) 

W^{u\z) 
S^{u\z) 



4i; 



^coth(2.-2.„)i/„, yiu\z) = j:-^^2u-2z.) 



a=l 

N 

^ sinh(2u - 2za) ' 

N 



N 



U^{u\z) = J2 



^±2{u-Za) 



sinh(2u - 20a) " 

a=l ^ 



V- 

■'-^ SIT 



25? 



sinh(2u — 2za) 
The corresponding semiclassical r-matrix has the form 

^7- I (g) r + cosh2u H®H + 2(5+ 5" + 5" 8) 5+) 
ih2u [3 



r{u) = 



sinhx 



+ e-2" [7- O C/+ + e2''C/+ (g) J/" + + VF" (g) W+] 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



Here we notice that (3.21) is equivalent to the r-matrix constructed out of the quadratic Casimir (3.14) 
in a standard way [18]. 

Substituting (3.21) and (3.19) into (2.4), we have 



Vl{u\z) ® l{v\z) + Vr{u - v) [l{u\z) (g) 1 + 1 (g) l{v\z)] 
= l{v\z) (g) l{u\z)P + [l{v\z) (g) 1 + 1 (g) l{u\z)] Vr{u - v), 



(3.22) 



whose consistence is guaranteed by the classical YB equation (1.3). 

^Prom (3.22) we can derive the commutation relations between the matrix elements of l{u\z). This 
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(2) 

gives us the defining relations of the A2 Gaudin algebra : 



[n{u\z),n{v\z 
[n{u\z),S'<'{v\z 

[n{u\z),W^{v\z 

[n{u\z),U'^{v\z), 

[y{u\z),yiv\z 

[y{u\z),W^{v\z 



[yiu\z 

[S^iu\z 
[S^{u\z 

[S^{u\z 

[S^{u\z) 
[U^{u\z 
[U^{u\z 

[U^{u\z) 

[U^iulz) 

[W^{u\z) 



,U^{v\z 
,S^{v\z 
,S^{v\z 

,U"^{v\z 

W^ivlz 
,U'^{v\z 
,U'^{v\z 

W'^ivlz 

W^{v\z 

W'^ivlz 



±2 

± 
± 



0, [niu\z),y{v\z)] = o, 

- cosh(2u - 2v)S'^{v\z) 
sinh(2ii — 2v) ' 
>VT(m|-z) - cosh(2u - 2v)W'^iv\z) 
sinh(2u — 2v) ' 
eT2(«-i;)^T(y,|2) - cosh(2M - 2v)U'^{v\z) 
sinh(2u — 2v) ' 
0, [y{u\z),S^iv\z)]=0, 
U^{u\z)-U^{v\z) 
sinh(2u — 2v) ' 

sinh(2u — 2v) 
0, [S^{u\z),W^iv\z)]=0, [S^{u\z),U^{v\z)] 
H{u\z) - n{v\z) 



T4 
±2 



sinh(2u — 2v) ' 

eT2(u-^)yy±(^|^) _ y\;±(v\z) 
sinh(2w — 2v) ' 
W^iulz) -W^ivlz) 



sinh(2u — 2v) ' 

0, [yV^{u\z),W'^{v\z)]=0, 
(H{u\z)-n{v\z)) 



^=F2(tt-«) 



T3 



sinh(2M — 2v) 

5^(M|z)-eT2("-^)5±(^;|z) 

sinh(2u — 2v) 
y{u\z) - e'^^^''-'"'^y{v\z) 



sinh(2w — 2v) 
n{u\z)-n{v\z) 



sinh(2u — 2v) 

A direct consequence of these relations is the commutativity of t^^\u\z) 

[t^'^\u\z),t^'^\v\z)]=0, yu,v 
from which the commutativity of the Gaudin Hamiltonians Ga follows immediately. 
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4 Off-shell Gaudin Equation 



In order to get the semiclassical limit of the OSBAE (2.14) we first consider the semiclassical expansions 
of the Bethe vectors defined in (2.13), (2.16) and (2.18): 

*„(ui,...,zi„|z) = (2r;)"$„(wi,...,M„|z) + o(7?"+i), 



^(n-i) - A2??j sinh(2n-2u,) '^n-i{u,\z)+o[r) ), 

*;i-2) = (2r/)"-^5-(n|^)$„_2(u;,n,» + o«), (4.1) 



where 



^n{ui,...,Un\z) = [W {Ui\z)+U {ui\z)] ^n-l{U2, ■■■,Un\z) 

g cosh(u,-t.,) ^""'^"'' ^^-'^ 

with $0 = |0) and ^i{ui\z) = [W-{ui\z) +U-{ui\z)] $o- 

The corresponding expansions of the c-number functions presented in the OSBAE are: (2.14) 



^j"-^) = 2r?/j"-^)+o(r?2), (4.4) 
^ cosh(itj - Mj) [ sinh(2u - 2ui) sinh(2u - 2uj) J ^ ' 

(4.5) 

where 

^ 2 
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A (2) AT ^ ^ ^ ^ ^ 

" " + 2" ~ 2 ^ COSh(w - Za)2 + 2 ^ COSh(u - U,)2 

JV n 

— {coth(M — Za) coth(M — Uj) + tanli(M — Za) tanli(M — Uj)} 

a=l j=l 

N 

+ {coth(M — Za) coth(u — Zb) + tanh(w — Za) tanh(u — Zb)} 

a<b 
n 

+ {coth(u — Uj) coth('u — Ufe) + tanh(u — Uj) tanh(u — Uk) 

j<k 

4 

+ - 



and 



(n-l) 



sinh(2u — 2uj) sinh(2w — 2uk) 



N 



(4.7) 



coth(Mj — Za) — {2 coth(uj — Uk) — tanh.{uj — Uk)} . (4.8) 

a=l kj^j 

Substituting these expressions into the Eq.(2.14) and comparing the coefficients of the terms 2(2j7)"+2 
we get the first non-trivial consequence for the semiclassical hmit of the OSBAE: 



n 2f^"~"^^6-' 

T^'\U\Z) $„(«!, Ur,\z) = A(2) ^^{u„ Ur,\z) - ^ ^""^^ 

■f^ smh(2u — 2uj) 



(4.9) 



Note that in this hmit the contributions from *(„_i) and ^'(4-2) combined to give a new vector 



valued function 



Z^-(M|z)e-2("-«^) + W-{u\z)] $„_i( uj\z) 



+S-{u\z) Yl 



k=l, k^j 



cosh(t(fe — Uj) 



$„_2( Uj,Uk\z), 



(4.10) 



Finally, we take the residue of (4.9) at the point u = Za to get the off-shell Gaudin equation: 

Ga^n{Ul,...,Un\z) = ga^„(Mi, M„|^) + V .'' ^^^T^, 

^ smh(2'Ui - 2za) 

a = 1,2, ...,N 



(4.11) 



where Qa is the residue of A 



(2) 



JV 



Qa = VeSu=Za-^n^ = ^ COt]l{Za - Zb) - COth(^;„ - Ul), 



1=1 



and is the residue of 



(4.12) 



= ^(W-e'"^-'^" +W-)$„_i(^,|z)+5-^ 



cosh(uj — Uk) 



^n-2{uk,Uj\z). (4.13) 
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The equation (4.11) allows us solve one of the main problems of the Gaudin models, i.e., the determi- 
nation of the eigenvalues and eigenvectors of the commuting Hamiltonians Ga (3.14): ga is the eigenvalue 
of Ga with eigenfunction provided ui are solutions of the following equations = 0, i.e.: 

N n 

coth(uj — Za) = {2 coth('Uj — Uk) — tanh(uj — Uk)} , j = 1,2,..., n. (4-14) 

o=l ki^j 

Moreover, as we will see in the next section, the off-shell Gaudin equation (4.11) provides solutions for 
the differential equations known as KZ equations. 

5 Knizhnik-Zamolodchickov equation 

The KZ differential equation 

= G'i(^)*(^), (5.1) 

appeared first as a holonomic system of differential equations of conformal blocks in WZW models of 
conformal field theory. Here ^'(^) is a function with values in the tensor product Fi • • • (g) Vat of 
representations of a simple Lie algebra, k = k + g , where k is the central charge of the associated 
Kac-Moody algebra, and g is the dual Coxeter number of the simple Lie algebra. 

One of the remarkable properties of the KZ system is that the coefficient functions Gi{z) commute 
and that the form w = Gi{z)dzi is closed [13]: 



dG, dG. 



dzi dzj 



[Gi,Gj]=0. (5.2) 



Indeed, it was indicated in [13] that the equations (5.2) are not just a flatness condition for the form oj 
but that the KZ connection is actually a commutative family of connections. 

In this section we will identify G, with our A2 Gaudin Hamiltonians Ga derived in the previous 



section 



JV 



^« = g sinh(2.' - 2z,) 1'^°^^^^^° - ^ + ^ 

+2 (5+ ® + S- ® S+) + W- W+ + W+ ® W,- 

a = 1,2,..., N. (5.3) 

and show that the corresponding differential equations (5.1) can be solved via the off-shell Bethe Ansatz 
method. 
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Let us now define the vector- valued function ^/{zi, zn) through multiple contour integrals of the 
Bethe vectors (4.2) 

^{zi, zn) ~ j' " ' j' ''^{'^\^)'^n{u\z)du\...dun, (5.4) 

where X {u\z) = X [ui, ...,Un,zi, ...,zn) is a scalar function which in this stage is still undefined. 
We assume that ^{zi, ...,zn) is a solution of the equations 

K ^— ^ = Ga*(zi,...,zw), a=l,2,...,iV (5.5) 

OZa 

where Ga are the Gaudin Hamiltonians (3.14) and k; is a constant. 
Substituting (5.4) into (5.5) we have 

where we are using a compact notation § {o} du = § . . . § {o} d,ui ■ ■ ■ d,u„. 

Using the Gaudin algebra (3.23) one can derive the following non-trivial identity 

dz^ ^ dui \^sinh(2w; - 2za) J' ^ ' ' 

which allows us write (5.6) in the form 

i^l X{u\z) . J du. (5.8) 

duiy ^ ' ^sinh(2w; -2^„)y ^ ' 

It is evident that the last term of (5.8) vanishes, because the contours are closed. Moreover, if the scalar 
function ^"(^12;) satisfies the following differential equations 

dX{u\z) dX{u\z) An-i)^, . . 

K = gaX{u\z), K = f] 'X{u\z), (5.9) 

we are recovering the off-shell Gaudin equation (4.11) from the first term in (5.8). 

Taking into account the definition of the scalar functions /j" ^''(4.8) and Qa (4.12), we can see that 
the consistency condition of the system (5.9) is insured by the zero curvature conditions jdza = 

dga/duj. Moreover, the solution of (5.9) is easily obtained 

N n N n 

X{u\z) = [] sinh(z„ - Zb)!/-^ [] [cosh(w, - Ufe)'^" sinh(w, - Ufc)-^/^] [] [] sinh(z„ - M,)-l/^ 

a#6 j#fe a=lj=l (5.10) 

This function determines the monodromy of ^'(^i, ■■■,Zn) as solution of the trigonometric KZ equation 
(5.5) and these results are in agreement with the Schechtman-Varchenko construction for multiple contour 
integral as solutions of the KZ equation in an arbitrary simple Lie algebra [12]. 
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6 Conclusion 

In this paper the A2 19-vertex model was investigated generalizing previous rational vertex models 
results relating the Gaudin magnet models to the semiclassical off-shell algebraic Bethe Ansatz of these 
vertex models. 

Using the semiclassical limit of the transfer matrix of the vertex model we derived the trigonometric 
A2 Gaudin Hamiltonians. The reduction of the off-shell Gaudin equation to an eigenvalue equation 
gives us the exact spectra and eigenvectors for the corresponding Gaudin magnets. Data of the off-shell 
Gaudin equation were used to show that a hypergeometric type integral (5.4) solves the trigonometric 
KZ differential equation. 

Our results corroborate the method already used with success to construct solutions of trigonometric 
KZ equations [19, 20, 21, 22] and elliptic KZ-Bernard equations [23]. 
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